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Modes  in  a Nonneutral  Plasma  Column 
of  Finite  Length 

S.  Neil  Rasband  and  Ross  L.  Spencer 


Department  of  Physics  and  Astronomy,  Brigham  Young  University,  Provo,  Utah  84602 


Abstract.  A Galerkin,  finite-element,  nonuniform  mesh  computation  of  the  mode 
equation  for  waves  in  a non-neutral  plasma  of  finite  length  in  a Cold-Fluid  model  gives 
an  accurate  calculation  of  the  mode  eigcnfrequencies  and  eigenfunctions.  We  report 
on  studies  of  the  following:  (l)finite-length  Trivelpiece-Gould  modes  with  flat-top  and 
realistic  density  profiles,  (2)finite-length  diocotron  inodes  with  flat  density  profiles. 
We  compare  with  the  frequency  equation  of  Fine  and  Driscoll  [Phys  Plasmas  5,  601 
(1998)]. 


INTRODUCTION 

The  familiar  Cold-Fluid  drift  model  for  the  nonneutral  plasma  gives  inside  the 
plasma  the  mode  equation  for  the  perturbed  potential  [Ij. 
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FIGURE  1.  Region  of  Computation 
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The  computation  region  is  illustrated  in  Figure  1 with  0 < r < ?Vaii  and  0 < ;:  < 
2waii,  where  the  plasma  in  this  region  is  confined  to  the  region  with  the  crosses. 
Equation  (1)  can  be  written  in  the  form 
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Construct  a decomposition  of  the  region  of  interest  into  triangular  elements, 
where  the  plasma  boundary  is  approximated  by  edyes  of  the  triangles.  Figure  2 
shows  an  example  with  rwaii  = 3.81cm  and  3w,,n  = 30cm. 


Quarter  cylinder  mesh 


-z(cm)- 


FIGURE  2.  A triangulatioii  of  a plasma  equilibrium.  Tlic  region  occupied  by  tlie  plasma  is 
shaded.  Note  that  the  scales  for  the  vertical  and  horizontal  axes  are  not  the  same. 


Each  triangle  has  6 nodes  (3  mid-points  for  the  sides  and  3 vertices)  and  on  each 
node  I a parabolic  function: 


y)  — /3i  + + p3ll  + + (hxy  + pel/. 


This  function  is  defined  so  that  it  has  value  1 at  the  7th  node  and  0 at  all  other 
nodes  in  the  triangle.  Then  approximate  as  a sum  over  nodes: 

¥'^.r,y)  = Y^Ci'^i{x,y). 

/ 
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Boundary  nodes  have  the  Cj  determined  from  boundary  conditions  on  y). 

The  Galerkin  integration  of  Eq.(l)  multiplied  by  the  approximating  functions 
Tj  proceeds  numerically  by  doing  one  triangular  element  at  a time.  If  the  element 
is  outside  the  plasma,  then  e = 1,  otherwise  it  is  as  given  in  Eq.(2).  This  gives  a 
matrix  equation  for  the  C/. 

Y^AjiCf  = 0 

/ 

with  nonzero  values  of  Cj  only  for  certain  values  (eigenvalues)  of  lo.  In  practice  we 
set 


'^AjiC[  = l,  for  each  J (3) 

I 

and  look  for  oj  such  that  max(C/)  — > oo  or  so  that  1/ max(C/)  — )•  0. 

TRIVELPIECE-GOULD  (M=0)  MODES 

As  a first  example  we  present  the  results  for  a flat-top  density  profile  with  the 
plasma  edge  at  rpiasma  = 1.89cm,  ^plasma  = 17.71cm.  The  triangulation  for  this 
equilibrium  is  shown  in  Figure  2.  The  aspect  ratio  a = 9.37  and  rpiasma/j’waii  = 
0.496.  For  the  modes  that  can  be  compared  with  the  results  in  Table  II  of  Jennings, 
Spencer,  and  Hansen  [2]  the  agreement  is  excellent. 

Figure  3 shows  a scan  in  frequency  for  even  modes  in  ^ with  some  of  the  prominent 
modes  indicated. 


(x10~^)  ''S  1/maxomp 


FIGURE  3.  Scan  of  l/max(C/)  as  described  following  Equation  (3)  for  even  modes. 
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A similar  scan  in  frequency  for  the  same  equilibrium  but  odd  modes  gives  fre- 
quencies w/wpe  = 0.1079,0.3060,0.4630,0.5766  for  the  modes  (1,0),  (3,0),  (5,0), 
and  (7,0),  respectively.  Figure  4 shows  the  perturbed  potential  eigenfunctions  for 
some  of  these  modes. 


Eigenfunction  ♦(r.z)  for  (2,0)  Mode  Eigenfunction  '^(r.z)  for  (4.0)  Mode 


Eigenfunction  ♦(r,z)  for  (1.0)  Mode  Eigenfunction  ♦(r.z)  for  (3,0)  Mode 


FIGURE  4.  perturbed  potential  eigenfunctions  for  selected  modes 


Lastly,  for  m=0  modes  we  consider  briefly  the  effect  of  a radial  dependence  in 
the  density  profile.  We  compute  e<iuilibria  [3]  whose  midplanc  density  is 

n(r)/rio  = (1  - \/(1.5/i/)(r/;Vaii)^) exp(-(7-/r,)‘'). 

We  choose  r;  = 7’waii/2  and  compare  two  choices  for  u,  u = 5.0,40.0,  corresponding 
to  a more-or-less  average  monotonic  profile  and  a flat  profile,  respectively.  Figure  5 
shows  the  profiles  with  the  corresponding  (2,0)  mode  frequencies.  The  changes  in 
the  mode  frequencies  for  such  changes  in  density  profiles  are  on  the  order  of  lO'X . 
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FIGURE  5.  Density  profiles  and  frequencies  for  (2,0)  inodes. 
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THE  DIOCOTRON  (M=l)  MODE 


We  examined  the  finite-length  diocotron  mode  frequency  for  a number  of  differing 
equilibria  with  varying  plasma  radii.  All  equilibria  have  flat-top  density  profiles 
and  are  computed  in  a Malmberg  trap  with  radius  7’waii  = 3.81cm,  half-length 
2waii  = 30.0cm  and  magnetic  field  375  G.  Figure  6 shows  the  shift  in  frequency 
from  the  infinite  length  result  as  a function  of  the  plasma  radius.  This  figure  also 
compares  these  results  to  those  obtained  with  the  formula  of  Equation  (24)  in  Fine 
and  Driscoll  [4]. 


Diocotron  Frequency  Shift  vs  Plosmo  Rodius 


FIGURE  6.  Diocotron  frequency  shift  as  a function  of  plasma  radius. 

For  the  case  of  Cwaii/fpiasma  = 2.0  Figure  7 shows  the  shape  of  the  perturbed 
plasma  potential  in  the  quarter  cylinder  computation  region.  For  a fixed  ^ value 
less  than  the  plasma  half-length  it  is  evident  that  inside  the  plasma  the  perturbed 
potential  dependence  on  r is  almost  linear.  A detailed  examination  of  these  eigen- 
functions shows  that  inside  the  plasma  their  dependence  on  r is  like  ?'  -|-  ar^  -!-■•• 
and  in  2 like  1 -f  bz^  -)-...  with  a and  b small  for  long  plasmas.  This  is  consistent 
with  solutions  inside  the  plasma  that  go  like  a modified  Bessel  function  /j  in  r and 
cosh(fcz)  in  z,  with  k very  small  corresponding  to  a wavelength  much  longer  that 
the  length  of  the  confining  cylinder  [1]. 

This  curvature  in  r is  readily  seen  for  a pancake-like  equilibrium.  In  Figure  8 
we  show  the  edge  curve  for  a pancake-like  equilibrium  with  Cpiasma  = 2.524cm  and 
^plasma  = 0.138cm  for  an  aspect  ratio  of  a = 0.055.  This  Figure  also  includes  a plot 
of  the  scaled  perturbed  potential  as  a function  of  r for  a fixed  value  of  2 where  the 
curvature  is  readily  apparent.  The  m=l  diocotron  mode  for  this  equilibrium  has  a 
frequency  of  w/wpe  = 0.00902. 
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Potentiol  <t(r,2) 


FIGURE  7.  Perturbed  Potential  in  the  quarter  cylinder  for  a m=l  diorotron  mode 


Plosmo  edge  curve  for  pancake  equilibrium 


•(r)  ot  a fixed  volue  of  z 
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FIGURE  8.  Pancake-like  plasma  edge  and  perturbed  potential  as  a function  of  »•  for  a fixed 
value  of  z = 0.05cm 
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